demo daphnia instructions
- Load the R scripts in chunk #1
- Run the chunk with the function demo_daphnia(). Run it from the console and not from the notebook otherwise it will take longer will not be able to see the animation.
# chunk #1
library(ReacTran)
library(simecol)
source("../R/daphnia_deb_phyto.R")
source("../R/daphnia_deb_equations.R")
source("../R/daphnia_deb_lifeequations.R")
source("../R/daphnia_deb.R")
source("../R/daphnia_deb_ibm.R")
source("../R/deb_ibm.R")

time= 1 ref. Vol= 10 abu= 10

time= 2 ref. Vol= 20 abu= 9.2

time= 3 ref. Vol= 20 abu= 8.85

time= 4 ref. Vol= 20 abu= 8.1

time= 5 ref. Vol= 20 abu= 7.8

time= 6 ref. Vol= 2.083333 abu= 49.92

time= 7 ref. Vol= 4.166667 abu= 47.04

time= 8 ref. Vol= 4.166667 abu= 44.64

time= 9 ref. Vol= 4.166667 abu= 42

time= 10 ref. Vol= 0.6127451 abu= 168.096

time= 11 ref. Vol= 1.22549 abu= 153.408

time= 12 ref. Vol= 1.22549 abu= 144.432

time= 13 ref. Vol= 1.22549 abu= 136.272

time= 14 ref. Vol= 0.4589851 abu= 235.3018

time= 15 ref. Vol= 0.9179702 abu= 215.6933

time= 16 ref. Vol= 0.9179702 abu= 162.3146

time= 17 ref. Vol= 0.9179702 abu= 124.187

time= 18 ref. Vol= 0.9179702 abu= 108.936

time= 19 ref. Vol= 2.753911 abu= 70.8084

time= 20 ref. Vol= 2.753911 abu= 39.21696

time= 21 ref. Vol= 5.507821 abu= 26.14464

time= 22 ref. Vol= 11.01564 abu= 14.16168

time= 23 ref. Vol= 22.03128 abu= 8.98722

time= 24 ref. Vol= 22.03128 abu= 5.08368

time= 25 ref. Vol= 66.09385 abu= 2.85957

time= 26 ref. Vol= 66.09385 abu= 1.57352

time= 27 ref. Vol= 198.2816 abu= 0.93806

time= 28 ref. Vol= 198.2816 abu= 0.57494

time= 29 ref. Vol= 396.5631 abu= 0.39338

time= 30 ref. Vol= 793.1262 abu= 0.2219067

time= 31 ref. Vol= 793.1262 abu= 0.1386917

time= 32 ref. Vol= 1586.252 abu= 0.08573667

time= 33 ref. Vol= 3172.505 abu= 0.05043333

time= 34 ref. Vol= 3172.505 abu= 0.04854208

time= 35 ref. Vol= 3172.505 abu= 0.04665083

time= 36 ref. Vol= 3172.505 abu= 0.04539

time= 37 ref. Vol= 3172.505 abu= 0.04318354

time= 38 ref. Vol= 3172.505 abu= 0.04129229

time= 39 ref. Vol= 3172.505 abu= 0.03877062

time= 40 ref. Vol= 3172.505 abu= 0.03687937

time= 41 ref. Vol= 3172.505 abu= 0.03435771

time= 42 ref. Vol= 3172.505 abu= 0.03246646

time= 43 ref. Vol= 6345.01 abu= 0.03026

time= 44 ref. Vol= 6345.01 abu= 0.02821115

time= 45 ref. Vol= 6345.01 abu= 0.02679271

time= 46 ref. Vol= 6345.01 abu= 0.02553187

time= 47 ref. Vol= 369.7558 abu= 0.2731532

time= 48 ref. Vol= 739.5116 abu= 0.2569263

time= 49 ref. Vol= 739.5116 abu= 0.2474606

time= 50 ref. Vol= 739.5116 abu= 0.2366427

time= 51 ref. Vol= 331.6196 abu= 0.2834573

time= 52 ref. Vol= 663.2391 abu= 0.2593333

time= 53 ref. Vol= 663.2391 abu= 0.2457635

time= 54 ref. Vol= 663.2391 abu= 0.2352093

time= 55 ref. Vol= 292.1758 abu= 0.3627952

time= 56 ref. Vol= 292.1758 abu= 0.3525274

time= 57 ref. Vol= 584.3517 abu= 0.3080337

time= 58 ref. Vol= 584.3517 abu= 0.3028998

time= 59 ref. Vol= 584.3517 abu= 0.2874981

time= 60 ref. Vol= 62.23127 abu= 1.783669

time= 61 ref. Vol= 62.23127 abu= 1.671185

time= 62 ref. Vol= 62.23127 abu= 1.606909

time= 63 ref. Vol= 124.4625 abu= 1.494425

time= 64 ref. Vol= 27.47518 abu= 3.748839

time= 65 ref. Vol= 54.95035 abu= 3.566856

time= 66 ref. Vol= 54.95035 abu= 3.403072

time= 67 ref. Vol= 54.95035 abu= 3.239288

time= 68 ref. Vol= 19.55529 abu= 5.420529

time= 69 ref. Vol= 19.55529 abu= 5.164844

time= 70 ref. Vol= 39.11057 abu= 4.909159

time= 71 ref. Vol= 39.11057 abu= 4.602336

time= 72 ref. Vol= 15.21812 abu= 6.636825

time= 73 ref. Vol= 6.503471 abu= 14.91511

time= 74 ref. Vol= 13.00694 abu= 14.60759

time= 75 ref. Vol= 13.00694 abu= 14.14629

time= 76 ref. Vol= 5.993982 abu= 16.1829

time= 77 ref. Vol= 2.814076 abu= 31.98208

time= 78 ref. Vol= 5.628152 abu= 29.84994

time= 79 ref. Vol= 5.628152 abu= 28.78387

time= 80 ref. Vol= 5.628152 abu= 27.36244

time= 81 ref. Vol= 1.803895 abu= 65.41401

time= 82 ref. Vol= 1.803895 abu= 63.19659

time= 83 ref. Vol= 1.803895 abu= 59.31609

time= 84 ref. Vol= 1.803895 abu= 58.20738

time= 85 ref. Vol= 1.803895 abu= 95.34923

time= 86 ref. Vol= 0.5352804 abu= 177.4771

time= 87 ref. Vol= 1.070561 abu= 162.5316

time= 88 ref. Vol= 1.070561 abu= 154.1248

time= 89 ref. Vol= 1.070561 abu= 144.7839

time= 90 ref. Vol= 1.070561 abu= 132.6408

time= 91 ref. Vol= 2.141122 abu= 88.73854

time= 92 ref. Vol= 2.141122 abu= 64.91925

time= 93 ref. Vol= 4.282243 abu= 41.567

time= 94 ref. Vol= 4.282243 abu= 29.89088

time= 95 ref. Vol= 8.564486 abu= 19.14884

time= 96 ref. Vol= 8.564486 abu= 12.25993

time= 97 ref. Vol= 25.69346 abu= 7.472719

time= 98 ref. Vol= 25.69346 abu= 4.787211

time= 99 ref. Vol= 51.38692 abu= 2.88011

time= 100 ref. Vol= 102.7738 abu= 1.517896

time= 101 ref. Vol= 205.5477 abu= 0.9243598

time= 102 ref. Vol= 205.5477 abu= 0.5302906

time= 103 ref. Vol= 616.643 abu= 0.2773079

time= 104 ref. Vol= 616.643 abu= 0.1670334

time= 105 ref. Vol= 1233.286 abu= 0.15406

time= 106 ref. Vol= 1233.286 abu= 0.143519

time= 107 ref. Vol= 1233.286 abu= 0.1362214

time= 108 ref. Vol= 1233.286 abu= 0.1329781

time= 109 ref. Vol= 1233.286 abu= 0.1248697

time= 110 ref. Vol= 1233.286 abu= 0.1183829

time= 111 ref. Vol= 1233.286 abu= 0.1151396

time= 112 ref. Vol= 1233.286 abu= 0.1086528

time= 113 ref. Vol= 1233.286 abu= 0.09811187

time= 114 ref. Vol= 1233.286 abu= 0.09081429

time= 115 ref. Vol= 1233.286 abu= 0.08594924

time= 116 ref. Vol= 2466.572 abu= 0.08027335

time= 117 ref. Vol= 1209.104 abu= 0.09097646

time= 118 ref. Vol= 145.325 abu= 0.6261827

time= 119 ref. Vol= 290.65 abu= 0.6124205

time= 120 ref. Vol= 290.65 abu= 0.5780148

time= 121 ref. Vol= 290.65 abu= 0.5573714

time= 122 ref. Vol= 74.90979 abu= 1.428385

time= 123 ref. Vol= 74.90979 abu= 1.334939

time= 124 ref. Vol= 149.8196 abu= 1.148048

time= 125 ref. Vol= 149.8196 abu= 1.101325

time= 126 ref. Vol= 43.0516 abu= 2.438934

time= 127 ref. Vol= 43.0516 abu= 2.346022

time= 128 ref. Vol= 86.10321 abu= 2.229882

time= 129 ref. Vol= 86.10321 abu= 2.171812

time= 130 ref. Vol= 30.00112 abu= 3.233213

time= 131 ref. Vol= 10.38101 abu= 9.632974

time= 132 ref. Vol= 20.76202 abu= 9.054996

time= 133 ref. Vol= 20.76202 abu= 8.814172

time= 134 ref. Vol= 20.76202 abu= 8.332523

time= 135 ref. Vol= 5.832028 abu= 17.6611

time= 136 ref. Vol= 11.66406 abu= 16.97523

time= 137 ref. Vol= 11.66406 abu= 16.46083

time= 138 ref. Vol= 11.66406 abu= 15.60349

time= 139 ref. Vol= 4.21085 abu= 22.56076

time= 140 ref. Vol= 8.421701 abu= 21.37335

time= 141 ref. Vol= 8.421701 abu= 20.66091

time= 142 ref. Vol= 8.421701 abu= 19.71098

time= 143 ref. Vol= 2.089752 abu= 47.37404

time= 144 ref. Vol= 1.029435 abu= 90.34086

time= 145 ref. Vol= 2.058869 abu= 87.42664

time= 146 ref. Vol= 2.058869 abu= 82.5696

time= 147 ref. Vol= 0.8990695 abu= 119.0119

time= 148 ref. Vol= 0.8990695 abu= 189.0844

time= 149 ref. Vol= 0.8990695 abu= 180.1863

time= 150 ref. Vol= 0.8990695 abu= 173.5127

time= 151 ref. Vol= 0.8990695 abu= 161.2779

time= 152 ref. Vol= 0.8990695 abu= 170.176

time= 153 ref. Vol= 0.8990695 abu= 114.5629

time= 154 ref. Vol= 1.798139 abu= 87.86863

time= 155 ref. Vol= 1.798139 abu= 68.40406

time= 156 ref. Vol= 1.798139 abu= 58.39371

time= 157 ref. Vol= 5.394417 abu= 33.36783

time= 158 ref. Vol= 5.394417 abu= 19.64995

time= 159 ref. Vol= 16.18325 abu= 11.86412

time= 160 ref. Vol= 16.18325 abu= 7.476866

time= 161 ref. Vol= 48.54975 abu= 4.078291

time= 162 ref. Vol= 48.54975 abu= 2.451094

time= 163 ref. Vol= 97.0995 abu= 1.503612

time= 164 ref. Vol= 194.199 abu= 0.8341958

time= 165 ref. Vol= 194.199 abu= 0.5200851

time= 166 ref. Vol= 582.597 abu= 0.288364

time= 167 ref. Vol= 1165.194 abu= 0.1287339

time= 168 ref. Vol= 2330.388 abu= 0.07981503

time= 169 ref. Vol= 2330.388 abu= 0.06865809

time= 170 ref. Vol= 2330.388 abu= 0.0622214

time= 171 ref. Vol= 2330.388 abu= 0.05964672

time= 172 ref. Vol= 2330.388 abu= 0.05664293

time= 173 ref. Vol= 2330.388 abu= 0.05492648

time= 174 ref. Vol= 2330.388 abu= 0.05321002

time= 175 ref. Vol= 2330.388 abu= 0.04977712

time= 176 ref. Vol= 2330.388 abu= 0.0523518

time= 177 ref. Vol= 2330.388 abu= 0.04891889

time= 178 ref. Vol= 2330.388 abu= 0.04763155

time= 179 ref. Vol= 2330.388 abu= 0.04720244

time= 180 ref. Vol= 2330.388 abu= 0.06822898

time= 181 ref. Vol= 1049.724 abu= 0.09335784

time= 182 ref. Vol= 294.0404 abu= 0.350292

time= 183 ref. Vol= 588.0809 abu= 0.3298866

time= 184 ref. Vol= 268.5301 abu= 0.3649498

time= 185 ref. Vol= 537.0601 abu= 0.3426059

time= 186 ref. Vol= 144.7601 abu= 0.6009941

time= 187 ref. Vol= 289.5203 abu= 0.5802702

time= 188 ref. Vol= 289.5203 abu= 0.5630003

time= 189 ref. Vol= 289.5203 abu= 0.5422763

time= 190 ref. Vol= 121.6472 abu= 0.8878134

time= 191 ref. Vol= 121.6472 abu= 0.8384904

time= 192 ref. Vol= 243.2944 abu= 0.780947

time= 193 ref. Vol= 243.2944 abu= 0.7357343

time= 194 ref. Vol= 86.89084 abu= 1.242939

time= 195 ref. Vol= 36.6628 abu= 2.863938

time= 196 ref. Vol= 73.32561 abu= 2.673009

time= 197 ref. Vol= 73.32561 abu= 2.522993

time= 198 ref. Vol= 26.28158 abu= 3.5386

time= 199 ref. Vol= 26.28158 abu= 5.441073

time= 200 ref. Vol= 26.28158 abu= 5.288875


# chunk #2
demo_daphnia <- function() {
## load individual-level model (DEB and phyto, differential equations)
submodel <- daphnia_deb_phyto()
## test one outer step alone
m <- sim(submodel)
out(m)
### load individual-based model (IBM, discrete, population dynamics)
deb_ibm <- daphnia_deb_ibm()
initfunc(deb_ibm) <- function(obj) {
p <- as.list(parms(obj))
parms(obj)["WAM"] <- length2weight(p$SAM, p$l2w) # weight at maturity
parms(obj)["WON"] <- length2weight(p$SON, p$l2w) # weight of neonates
age <- sample(1:30, 100, repl=TRUE)
len <- p$SON + 2 * age/30 ## simple linear approximation
weight <- length2weight(len, p$l2w)
initlist <- list(X = c(x1=0, x2=0, x3=0.25),
Z = data.frame(
age = age,
weight = weight,
egg = sample(1:5, 100, repl=TRUE),
eggage = runif(100, min=0, max=5),
length = len,
eggstorage = 0
),
refVol = 10 # 10L
)
init(obj) <- initlist #list2vec(initlist)
obj
}
observer(deb_ibm) <- function (statelist) {
X <- statelist$X
Z <- statelist$Z
refVol <- statelist$refVol
hist(Z$length, col="red", breaks=seq(0.5,3,0.25), ylim=c(0,150))
c(X["x1"], X["x2"], X["x3"], abundance = nrow(Z)/refVol,
age = mean(Z$age), weight = mean(Z$weight), eggs = mean(Z$egg),
length = mean(Z$length), refVol = refVol)
}
#deb_ibm <- initialize(deb_ibm)
times(deb_ibm)["to"] <- 200
parms(deb_ibm)$samplesize <- c(100, 200)
parms(deb_ibm)$life.span <- 60
deb_ibm <- sim(deb_ibm)
o <- out(deb_ibm)
### plot results
par(mfrow=c(3,1))
plot(o$time, o$x1, type="l",
main="Phytoplankton", col="red", ylim=c(0,0.5),xlab="Day",ylab="mg C / L")
lines(o$time, o$x2, col="blue")
lines(o$time, o$x3, col="green")
plot(o$time, o$abundance, type="l", main="Abundance",xlab="Day", ylab="Ind / L")
bm.mgc <- o$weight * o$abund / 1000
plot(o$time, bm.mgc , type="l", main="Biomass",xlab="Day", ylab="mg C / L")
}
demo_daphnia()

time= 1 ref. Vol= 10 abu= 10

time= 2 ref. Vol= 20 abu= 9.8

time= 3 ref. Vol= 20 abu= 9.3

time= 4 ref. Vol= 20 abu= 8.8

time= 5 ref. Vol= 20 abu= 8.1

time= 6 ref. Vol= 2.053388 abu= 49.187

time= 7 ref. Vol= 4.106776 abu= 46.752

time= 8 ref. Vol= 1.740159 abu= 59.76464

time= 9 ref. Vol= 3.480319 abu= 56.31668

time= 10 ref. Vol= 0.6481041 abu= 160.4681

time= 11 ref. Vol= 0.6481041 abu= 154.2962

time= 12 ref. Vol= 0.6481041 abu= 171.2688

time= 13 ref. Vol= 0.6481041 abu= 160.4681

time= 14 ref. Vol= 0.6481041 abu= 197.4991

time= 15 ref. Vol= 0.6481041 abu= 185.1555

time= 16 ref. Vol= 1.296208 abu= 138.8666

time= 17 ref. Vol= 1.296208 abu= 101.8355

time= 18 ref. Vol= 2.592416 abu= 74.83366

time= 19 ref. Vol= 2.592416 abu= 46.6746

time= 20 ref. Vol= 5.184832 abu= 30.4735

time= 21 ref. Vol= 5.184832 abu= 21.02286

time= 22 ref. Vol= 15.5545 abu= 12.53657

time= 23 ref. Vol= 15.5545 abu= 7.39336

time= 24 ref. Vol= 31.10899 abu= 4.436016

time= 25 ref. Vol= 62.21799 abu= 2.507313

time= 26 ref. Vol= 124.436 abu= 1.52689

time= 27 ref. Vol= 124.436 abu= 0.892025

time= 28 ref. Vol= 373.3079 abu= 0.5303932

time= 29 ref. Vol= 373.3079 abu= 0.3402017

time= 30 ref. Vol= 746.6159 abu= 0.206264

time= 31 ref. Vol= 746.6159 abu= 0.1366164

time= 32 ref. Vol= 1493.232 abu= 0.1178652

time= 33 ref. Vol= 1493.232 abu= 0.1111683

time= 34 ref. Vol= 1493.232 abu= 0.1051411

time= 35 ref. Vol= 1493.232 abu= 0.0984442

time= 36 ref. Vol= 1493.232 abu= 0.09040794

time= 37 ref. Vol= 1493.232 abu= 0.08438074

time= 38 ref. Vol= 1493.232 abu= 0.08371105

time= 39 ref. Vol= 1493.232 abu= 0.07567479

time= 40 ref. Vol= 1493.232 abu= 0.06830822

time= 41 ref. Vol= 2986.464 abu= 0.0636204

time= 42 ref. Vol= 2986.464 abu= 0.06161133

time= 43 ref. Vol= 669.6107 abu= 0.1403801

time= 44 ref. Vol= 1339.221 abu= 0.1344065

time= 45 ref. Vol= 195.5068 abu= 0.5114913

time= 46 ref. Vol= 391.0135 abu= 0.4859167

time= 47 ref. Vol= 186.1969 abu= 0.5639191

time= 48 ref. Vol= 186.1969 abu= 0.5370658

time= 49 ref. Vol= 80.25729 abu= 1.183693

time= 50 ref. Vol= 160.5146 abu= 1.146313

time= 51 ref. Vol= 77.9197 abu= 1.360375

time= 52 ref. Vol= 77.9197 abu= 1.334707

time= 53 ref. Vol= 77.9197 abu= 1.732553

time= 54 ref. Vol= 77.9197 abu= 1.668384

time= 55 ref. Vol= 77.9197 abu= 1.540047

time= 56 ref. Vol= 77.9197 abu= 2.027729

time= 57 ref. Vol= 77.9197 abu= 2.41274

time= 58 ref. Vol= 27.92821 abu= 3.75964

time= 59 ref. Vol= 27.92821 abu= 3.652221

time= 60 ref. Vol= 27.92821 abu= 4.798016

time= 61 ref. Vol= 27.92821 abu= 4.511568

time= 62 ref. Vol= 10.01011 abu= 9.9899

time= 63 ref. Vol= 20.02022 abu= 9.790102

time= 64 ref. Vol= 8.483144 abu= 12.14172

time= 65 ref. Vol= 8.483144 abu= 11.78808

time= 66 ref. Vol= 8.483144 abu= 21.80795

time= 67 ref. Vol= 8.483144 abu= 20.62914

time= 68 ref. Vol= 8.483144 abu= 23.57616

time= 69 ref. Vol= 3.70443 abu= 28.0745

time= 70 ref. Vol= 3.70443 abu= 39.41227

time= 71 ref. Vol= 3.70443 abu= 41.3019

time= 72 ref. Vol= 3.70443 abu= 40.22211

time= 73 ref. Vol= 3.70443 abu= 48.86042

time= 74 ref. Vol= 1.722991 abu= 56.87785

time= 75 ref. Vol= 1.722991 abu= 71.96789

time= 76 ref. Vol= 1.722991 abu= 67.90519

time= 77 ref. Vol= 1.722991 abu= 88.2187

time= 78 ref. Vol= 1.722991 abu= 84.73639

time= 79 ref. Vol= 0.640517 abu= 131.1441

time= 80 ref. Vol= 1.281034 abu= 129.5828

time= 81 ref. Vol= 1.281034 abu= 150.6596

time= 82 ref. Vol= 1.281034 abu= 124.8991

time= 83 ref. Vol= 1.281034 abu= 107.7255

time= 84 ref. Vol= 2.562068 abu= 67.13327

time= 85 ref. Vol= 2.562068 abu= 46.44686

time= 86 ref. Vol= 5.124136 abu= 30.44416

time= 87 ref. Vol= 5.124136 abu= 21.07673

time= 88 ref. Vol= 15.37241 abu= 10.14805

time= 89 ref. Vol= 30.74481 abu= 5.854646

time= 90 ref. Vol= 30.74481 abu= 3.480262

time= 91 ref. Vol= 61.48963 abu= 2.276807

time= 92 ref. Vol= 122.9793 abu= 1.39861

time= 93 ref. Vol= 122.9793 abu= 0.8375396

time= 94 ref. Vol= 368.9378 abu= 0.4960186

time= 95 ref. Vol= 368.9378 abu= 0.3903097

time= 96 ref. Vol= 368.9378 abu= 0.3659154

time= 97 ref. Vol= 368.9378 abu= 0.3577839

time= 98 ref. Vol= 368.9378 abu= 0.3442315

time= 99 ref. Vol= 368.9378 abu= 0.3333896

time= 100 ref. Vol= 368.9378 abu= 0.3062847

time= 101 ref. Vol= 368.9378 abu= 0.2981533

time= 102 ref. Vol= 368.9378 abu= 0.2981533

time= 103 ref. Vol= 368.9378 abu= 0.2873113

time= 104 ref. Vol= 368.9378 abu= 0.2710484

time= 105 ref. Vol= 368.9378 abu= 0.3875992

time= 106 ref. Vol= 368.9378 abu= 0.4445194

time= 107 ref. Vol= 368.9378 abu= 0.4499404

time= 108 ref. Vol= 368.9378 abu= 0.425546

time= 109 ref. Vol= 41.54705 abu= 2.238426

time= 110 ref. Vol= 83.09409 abu= 2.118081

time= 111 ref. Vol= 40.5337 abu= 2.368399

time= 112 ref. Vol= 81.06741 abu= 2.09702

time= 113 ref. Vol= 25.57331 abu= 4.692391

time= 114 ref. Vol= 25.57331 abu= 4.457772

time= 115 ref. Vol= 25.57331 abu= 5.083424

time= 116 ref. Vol= 25.57331 abu= 4.770598

time= 117 ref. Vol= 25.57331 abu= 7.625136

time= 118 ref. Vol= 25.57331 abu= 7.507826

time= 119 ref. Vol= 11.16739 abu= 9.491926

time= 120 ref. Vol= 11.16739 abu= 9.13374

time= 121 ref. Vol= 22.33477 abu= 8.417368

time= 122 ref. Vol= 4.370797 abu= 24.70945

time= 123 ref. Vol= 4.370797 abu= 24.70945

time= 124 ref. Vol= 4.370797 abu= 25.85341

time= 125 ref. Vol= 4.370797 abu= 24.93824

time= 126 ref. Vol= 1.783999 abu= 56.61439

time= 127 ref. Vol= 3.567998 abu= 53.8117

time= 128 ref. Vol= 1.585777 abu= 61.79937

time= 129 ref. Vol= 3.171553 abu= 61.16876

time= 130 ref. Vol= 1.026393 abu= 110.0943

time= 131 ref. Vol= 1.026393 abu= 100.3515

time= 132 ref. Vol= 2.052785 abu= 96.45431

time= 133 ref. Vol= 2.052785 abu= 91.09574

time= 134 ref. Vol= 0.986916 abu= 95.2462

time= 135 ref. Vol= 0.986916 abu= 132.7367

time= 136 ref. Vol= 0.986916 abu= 112.4716

time= 137 ref. Vol= 1.973832 abu= 85.11362

time= 138 ref. Vol= 1.973832 abu= 70.42139

time= 139 ref. Vol= 1.973832 abu= 54.20927

time= 140 ref. Vol= 5.921496 abu= 27.86458

time= 141 ref. Vol= 5.921496 abu= 17.0565

time= 142 ref. Vol= 11.84299 abu= 11.31471

time= 143 ref. Vol= 23.68598 abu= 6.332859

time= 144 ref. Vol= 47.37197 abu= 4.137468

time= 145 ref. Vol= 47.37197 abu= 2.70202

time= 146 ref. Vol= 94.74394 abu= 1.519886

time= 147 ref. Vol= 189.4879 abu= 0.9499288

time= 148 ref. Vol= 189.4879 abu= 0.8813229

time= 149 ref. Vol= 189.4879 abu= 0.8602133

time= 150 ref. Vol= 189.4879 abu= 0.8391038

time= 151 ref. Vol= 189.4879 abu= 0.7599431

time= 152 ref. Vol= 189.4879 abu= 0.717724

time= 153 ref. Vol= 189.4879 abu= 0.6649502

time= 154 ref. Vol= 189.4879 abu= 0.6385633

time= 155 ref. Vol= 189.4879 abu= 0.6121764

time= 156 ref. Vol= 189.4879 abu= 0.5857895

time= 157 ref. Vol= 189.4879 abu= 0.5277382

time= 158 ref. Vol= 189.4879 abu= 0.8496586

time= 159 ref. Vol= 189.4879 abu= 0.8179943

time= 160 ref. Vol= 189.4879 abu= 0.7704978

time= 161 ref. Vol= 21.03084 abu= 3.899036

time= 162 ref. Vol= 42.06168 abu= 3.61374

time= 163 ref. Vol= 42.06168 abu= 3.471093

time= 164 ref. Vol= 42.06168 abu= 3.233347

time= 165 ref. Vol= 13.48131 abu= 7.269324

time= 166 ref. Vol= 26.96262 abu= 6.750087

time= 167 ref. Vol= 26.96262 abu= 6.342114

time= 168 ref. Vol= 26.96262 abu= 6.156672

time= 169 ref. Vol= 10.29108 abu= 9.717158

time= 170 ref. Vol= 20.58215 abu= 9.134128

time= 171 ref. Vol= 20.58215 abu= 9.5714

time= 172 ref. Vol= 20.58215 abu= 9.085543

time= 173 ref. Vol= 7.73765 abu= 15.63782

time= 174 ref. Vol= 2.803497 abu= 36.02644

time= 175 ref. Vol= 5.606993 abu= 34.59965

time= 176 ref. Vol= 5.606993 abu= 32.99451

time= 177 ref. Vol= 5.606993 abu= 31.03268

time= 178 ref. Vol= 1.704253 abu= 62.7841

time= 179 ref. Vol= 3.408506 abu= 57.50319

time= 180 ref. Vol= 3.408506 abu= 54.86274

time= 181 ref. Vol= 3.408506 abu= 53.10244

time= 182 ref. Vol= 0.8876319 abu= 112.6593

time= 183 ref. Vol= 1.775264 abu= 108.1529

time= 184 ref. Vol= 0.6699109 abu= 155.2445

time= 185 ref. Vol= 1.339822 abu= 146.2881

time= 186 ref. Vol= 0.4481009 abu= 229.8589

time= 187 ref. Vol= 0.4481009 abu= 263.3335

time= 188 ref. Vol= 0.8962018 abu= 218.7007

time= 189 ref. Vol= 0.8962018 abu= 181.8787

time= 190 ref. Vol= 0.8962018 abu= 152.8674

time= 191 ref. Vol= 1.792404 abu= 101.5396

time= 192 ref. Vol= 1.792404 abu= 56.34892

time= 193 ref. Vol= 5.377211 abu= 30.12714

time= 194 ref. Vol= 5.377211 abu= 19.15491

time= 195 ref. Vol= 16.13163 abu= 10.78626

time= 196 ref. Vol= 16.13163 abu= 7.438801

time= 197 ref. Vol= 32.26327 abu= 4.21532

time= 198 ref. Vol= 64.52653 abu= 2.820545

time= 199 ref. Vol= 64.52653 abu= 1.54975

time= 200 ref. Vol= 193.5796 abu= 0.9453476

